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ABSTRACT 



Gaia mission will offer an exceptional opportunity to per- 
form variability studies. The data homogeneity, its op- 
timised photometric systems, composed of 11 medium 
and 4-5 broad bands, the high photometric precision in G 
band of one milli-mag for V = 13—15, the radial veloc- 
ity measurements and the exquisite astrometric precision 
for one billion stars will permit a detailed description of 
variable objects like stars, quasars and asteroids. 

However the time sampling and the total number of mea- 
surements change from one object to an other because of 
the satellite scanning law. The data analysis is a challenge 
because of the huge amount of data, the complexity of the 
observed objects and the peculiarities of the satellite, and 
needs a thorough preparation. 

Experience can be gained by the study of past and present 
survey analysis and results and Gaia should be put in 
perspective with the future large scale surveys, like Pan- 
STARRS or LSST 

We present the activities of the Variable Star Working 
Group and a general plan to digest this unprecedented 
data set, focusing here on the photometry. 
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1. INTRODUCTION 



The logarithmic representation of the Universe by 
LGottetal.'('200l) 

(http://www.astro.princeton.edU/~mjuric/universe/l 
clearly shows the breadth and diversity of the Gaia 
mission Science. 

Gaia will measure targets from the solar system small 
bodies to the distant quasars, the bulk of objects being 
1 billion stars in our Galaxy. 

A huge number of the objects will be variable. Indeed 
asteroids are generally variable at the level SV ^ 0. 1-0.7 
mag. Though the OSOs are sometimes less variable than 
usually thought jSchecht er 2005), a very large fraction of 
them will be detected as variable with Gaia. 



On the side of stars, there is a Zoo of variables: If we 
look at the GCVS we may count more that 100 types and 
subtypes. It is customary to divide the variable stars in 
five main classes: Pulsating stars, eclipsing binaries, vari- 
ability induced by rotation, eruptive stars and cataclysmic 
variable stars. The time scale and amplitude of the vari- 
ation may be very different: from short event of seconds 
to secular evolution and from milli-magnitudes to several 
magnitudes. The variability may be periodic, or irregular 

The several instrumental devices on board, photometric, 
astrometric and radial velocity, will allow a deep investi- 
gation on variability induced phenomena. Gaia will allow 
a global and precise description of stellar variability in 
terms of physical quantities (temperature, gravity, metal- 
licity). On classical topics, Gaia will give access to abso- 
lute masses and si zes for about 1 000 eclipsing binaries 
as pointed out by Zwitter ( 2003), the period-luminostiy- 
metallicity realtions of Cepheids will greatly benefit from 
the parallaxes (cf. Ngeow & Kanbur, these proceedings). 
Asteroseismology will also benefit from accurate paral- 
laxes, because luminos ity and radius are stringent test for 
the stellar structure, cf. iFavatal (Il999l) . 

The G-band photometric and astrometric informations 
are deduced from the same CCD counts, therefore an in- 
terdependence between the two processings exists. As 
observed in other surveys, variability may affect astro- 
metric measurements (cf. £ourbaix et al. 2003) and 
astrometry (proper motions) may affect phot ometry by 
the reduction methods, cf. lEver & WozniakI (2001) or 
[Alcock et al. (2001b). Non uniform light distribution in 
Long Period Variables (LPVs) may affect the measured 
astrometry through photocenter wobbling (Jorissen, pri- 
vate communication). 

For these reasons, the variability analysis is an unavoid- 
able step of the Gaia mission to improve both scientific 
returns and astrometric quality. 



2. PAST SURVEYS 

In recent years, only the HIPPARCOS satellite has done 
a multi-epoch photometric all-sky survey with an associ- 
ated analysis aimed at systematically detecting variabil- 
ity. The HIPPARCOS programme consists of 1 18 204 stars 
among them a survey of 52 045 stars, complete down to 
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a well defined limiting magnitude function of the star 
colour and its ecliptic latitude. Results on variable stars 
were published in the volum es 1 1 and 12 and among the 
flags of the main catalogue (ESA 1997). The rate of vari- 
able stars is very high, about 10% with a similar rate in 
the survey and in the preselected star sample. 

Microlensing surveys like EROS, MACHO, OGLE de- 
tected and characterised a large number of variable stars 
with a lesser detection rate of about 1%. To mention one, 
OGLE-II observed about 40 million stars jUdalski et ail 
(1998), Udalskietal. (2000), lUdalski eUl (2002)) and 
IZebrunetal. (2001) Wozniak et al. (2002) extracted 
about 290 000 variable stars, data publicly available. 

From OGLE and HIPPARCOS surveys, we can see that 
there is a certain time lag between the data acquisition 
end and the published results of the variability. The 
OGLE-II survey collected data from 1997 to 2000, pub- 
lished the photometry for the variable objects in 2001 
for the Magellanic clouds and in 2002 for the bulge. In 
2004, the database is not fully scientifically exploited yet. 
HiPPARCOS finished the measurements in 1993, the cat- 
alogue was published in 1997, though the full (prepared) 
variability analysis was completed within a few months. 

These examples show that a thorough preparation of the 
analysis should be made well in advance in order not to 
lag with the data flow. 



3. ONGOING AND FUTURE SURVEYS 



The study of variability for Gaia should also be 
put in perspective with other projects. A website 
lists about 90 ex isting or future robotic telescopes 
( |http : // alpha.uni- sw .gwdg.de/ ~hessman/MONET/ 1 . 

We mention here three ongoing surveys: 

- ASAS (All-Sky Automated Survey, lPoimanski|[T997h 
is a project which goal is the photometric monitoring of 
approximately 10 million stars brighter than magnitude 
14. The photometry is in V and / bands. The data for 
variable stars is publicly available. There is an alert sys- 
tem. 

- OGLE (Optical Gravitational Lensing Experiment) 
is continuing to observe the bulge and the Magellanic 
clouds. It has also dedicated campaigns aiming at the de- 
tection of planetary transits (Udalski et al. 2003). OGLE 
is the first team to have detected planetary systems with 
the tran s it met hod. OGLE-III has a real time analysis 
lirdalskilil2(K)l . 

- SDSS (Sloan Digital Sky Survey) provides multi-epoch 
measurements. In the Northern sky, the number of epochs 
is limited to two but thanks to the multi-colour data, 
it allows to detect efficiently some va riability types cf. 
Ilvezic et alJ (l2000h . Ilvezic etal.1 ll2003h . We note that in 
the Southern Survey the number of epochs will be about 
a dozen. 

There are two large scale future surveys which should be 
mentioned: 



- Pan-STARRS (Panoramic Survey Telescope & Rapid 
Response System) will monitor the sky down to a limiting 
magnitude of 24 about 20 times per year. The survey will 
be executed by four 1.8m telescopes located in Hawaii. 
Each telescope will have a 3 deg field of view. The first 
light of the first telescope is planned for 2006 and the full 
system will be operational in 2008. The photometric sys- 
tem will have the Sloan g, r, i bands, a "sloanish" z band 
and a y band at 1 .7 micron. The photometric precision at 
the bright end will be about 5 mmag. 

- LSST (Large Synoptic Survey Telescope) will cover 
down to a limiting magnitude of 24 in 5 bands 20000 
square degrees with a single 8.4m telescope. The first 
light is planned for 2012. 



4. PHOTOMETRIC SURVEYS FROM SPACE 



Two space missions, COROT and Kepler, will describe 
as by products the behaviour of large numbers of small 
amplitude variable stars thanks to their high photomet- 
ric precision. Two remar kable small sa tellites, WIRE 
dluzasi 2002) and MOST ( Walker et al. 2003), have been 
exploring the small amplitude variability domain working 
both in a single target mode: 

- The CNES/ESA COROT mission (http:// 
smsc.cnes.fr/COROT/index.htm): The satellite will 
be in operation from 2006 to 2009 with a possible 
extension. There will be two programs run in parallel, 
one dedicated to asteroseismology, the other to planetary 
transit searches. The exo-planet program will monitor 
in total 30000-60000 stars down to V=15.5 with colour 
information, in five different fields, during 150 days 
each. Measurements will be transferred to Earth every 
8 minutes. After 1 hour of integration, mean intensities 
will have a precision of 700 ppm (parts per million) 
at V — 15.5. The asteroseismology program will 
monitor, 5 different fields containing one primary target 
per field brighter than V — 6 and a maximum of 9 
secondary targets during 150 days. Finally there will be 
4 exploratory programs of 30 days each. 

The Kepler Mission (NASA, http:// 
www.KeplerNASA.gov/): It will be launched in 
2007 and will be operational for 4 years, with a possible 
2 year mission extension. At least 100000 preselected 
main sequence stars, in the magnitude range from 9 
to 15, will be continuously monitored thanks to the 
Earth-trailing heliocentric orbit and the position of the 
observed field selected to avoid interruptions in the 
data acquisition. After 6.5 hours of integration, the 
mean intensity, for a source of mag V=12, will have a 
precision of <14 ppm (including instrument, background 
and photon shot noises). Additionally, 200 sources 
win be monitored for asteroseismological purposes at 
a one-minute cadence. These source will be changed 
every 1 to 3 months so as to have monitored all the bright 
dwarfs in the FOV for p- and g-modes, about 5 000 of 
them. 
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5. GAIA SURVEY 



Gaia will observe 1 billion objects with a mean of 
^80 measurements in the Astro-field (G-band and Broad 
Band Photometry-BBP), and ^90 measurements in the 
Medium Band Photometry (MBP) during 5 years. Possi- 
bly a near IR magnitude from the Spectro-field can also 
be computed. The photometric systems are still under 
discussions. The photometric precision in the G-band is 
very high from 1 mmag at V=13-15 to 0.02 mag at the 
limiting magnitude V=20; for more details on the photo- 
metric performances, see Jordi ( 2004) . 



5.1. Gaia Sampling 

The satellite will be operated at the second Lagrangian 
point. The sampling law is governed by two main factors: 

1. The rotation of the satellite on itself in 6 hours. 

2. The precession of the satellite rotation axis on a cone 
of 50 deg angle to the Sun, in 70 days. 

The Astro-field will gather two observing directions sep- 
arated by an angle of 106 degrees: the preceding (PFOV) 
and following (FFOV) fields of view. In the Astro-field, 
a star has measurements generally grouped consecutively 
with the time intervals of lh46, 4hl4, lh46, etc, corre- 
sponding to observations in the PFOV-FFOV-PFOV, etc. 
Then the next group will occur typically 30 days later 
Both the total length of the sequence of short consecutive 
measurements and the time intreval between the grouped 
measurements depend on the position of the star on the 
celestial sphere. We show in Figure^the dependence on 
the number of measurements as function of the ecliptic 
latitude /?. The mean number of measurements is 82 but 
varies between 45 to 210. 
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Figure 1. Number of measurements (transits) as function 
of the ecliptic latitude [3. The mean number of measure- 
ments is 82. 



5.2. Number of variable objects detected by Gaia 

The number of variable objects detected by Gaia 
is still quite uncertai n. First estimations done by 
lEver & Cuvp"ersl ll2000h lead to about 18 million classi- 
cal variable stars. 

Caution is still needed: for instance the estimates of the 
number of eclipsing binaries by different authors differ 
greatly: Soederhjelm (private communication) estimated 
this number to be half a million, whereas Ever & Cuvper^ 
(12000) estimated it to 2-3 milHons and Zwitter ( 2002) to 7 
millions.We list several estimates for other variable types: 

- lEver&CuvpersI MmCf) estimated the numbers of RR 
Lyr to 70 000, of Cepheids to 2 000-8 000, of Mira-LPVs 
to 200000, of 5 Scuti stars to 60000-240000, of SPBs to 
15 000, of (3 Cep stars to 3 000. 

- Robichon (2002, private communication) estimated the 
number of planetary transit systems to 5 000-30 000. 
-iBelokurov & Evans (2003) estimated the number of Su- 
p ernovae to 2 1 000. 

- iTotani & Panaitescul(l2002h estimated that Gaia will ob- 
serve 720 optical counter parts of gamma ray bursts. 



5.3. Rate of correct detection of Periods 

Eyer and Mignard has investigated the rate of correct de- 
tection of the period for a given periodic signal observed 
with the Gaia sampling law. Sinusoidal signals have 
been generated for pairs of P (period) and Q (amplitude 
to noise ratio) and sampled with the Gaia scanning law 
for 10000 stars uniformly distributed over the celestial 
sphere. For each star a periodogram is constructed, and 
the highest power period Ph is confronted to the true pe- 
riod P. About 1 .4 million periodograms have been com- 
puted. 

As expected, the rate of correct period detection depends 
on the number of measurements and therefore on the 
ecliptic latitude. Once the number of measurements ef- 
fect is corrected, the dependence on the ecliptic latitude 
does not fully vanish. A dip reaching 10-15% in the pe- 
riod correct detection rate is still present at ecliptic lati- 
tude 20 degrees. 

Very high amplitude to noise ratio Q is not necessary to 
recover the period. Already at Q ~ 1.3, the rate of correct 
detection is mostly above 90%. For example a sinusoidal 
signal of amplitude 3 mmag ai V = 14, will have its 
period correctly recovered from several hours to hundreds 
of days. 



5.4. Variable Star Working Group activities 

In the previous sections we covered some activities of the 
Working Group. They are fully described at the following 
website: 

\protect\vrule widthOpt \protect \href { http : //obswww 
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You may find more information on: 

- Variable star characteristics: We prepare a uniform de- 
scription of variable star types. Information collected are: 
general references, physical properties, variability prop- 
erties (times scales and amplitudes of the photometric and 
radial velocity variations), physical source of variability, 
number of star known, fraction of variable is the same 
area of the HR diagram. Once the filters of Gaia will be 
defined, variability at the filter wavelengths will be esti- 
mated for each type. Such characteristics will serve for 
completing the Galactic Model used for GDAAS (Gaia 
Data Access and Analysis Study), for the estimation of 
the number of variables, as well as for the classification. 

- Compilations of past, present and future surveys: We 
want to keep track of the research done on these surveys. 
Therefore we give a list of past surveys with their returns 
on variable stars. We will add the methods used to anal- 
yse the data. We will also give the characteristics of some 
future surveys. 

- Forecasts about variable stars to be detected by 
Gaia: We gather published information from different re- 
searchers on the forecasts. 

- Tools: We give a list of typical tools which are used in 
variable star research. Algorithms of Period search, time 
scales, or classification methods. We put online software 
for some algorithms, or links to the codes whenever pos- 
sible. 

- GDAAS related activities: We are developing software 
for the GDAAS analysis (see the following section). 



6. GLOBAL ANALYSIS 



The variability analysis can be seen as a process of data 
compression. If a star has a constant magnitude, then the 
whole time series can be summarised by two quantities 
(per band): the mean magnitude and its precision. For 
instance the compression factor for the G-band photome- 
try with 80 measurements is about 140 (the measurement 
epochs can be removed). If the object is detected as vari- 
able, then the task to compress the data is more labori- 
ous, the time series may be reduced, for example, to am- 
plitude(s), time-scale(s), period(s), epochs, parameters of 
Fourier series, variable types. The goals of the global 
analysis are to separate the observed objects in constant 
(or more precisely non detected as variable) and variable, 
and to reach a classification with the associated relevant 
parameters. We proposed an organigram to fulfil these 
goals in Figure|2l 



6.1. Detection of variability 

The first step is to detect variability. As different vari- 
ability behaviours have different photometric signatures, 
a set of different statistical tests had to be designed. 

Thanks to the very high G-magnitude precision, it has 
been decided, during the VSWG Cambridge meeting in 
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Figure 2. Organigram of used data processing of the 
variable stars, baseline of the current development of the 
software for GDAAS. 



2004, to use the magnitudes deduced from one passage of 
the observed object over the 1 1 CCDs of the Astro-field. 
Studies of very short term variability should also be ex- 
plored using the per CCD information (most algorithms 
developed for the per transit analysis are applicable at this 
level). 

To detect variability we perform a number of statistical 
tests on the G-band photometric data. For these tests two 
aspects are important: 

1 . they should be mathematically formulated in a cu- 
mulative way, so that during the mission, when new 
data are coming, the values of the statistical tests can 
be updated from old cumulated quantities and the 
new data. So, the variability status of every object 
can be monitored. 

2. for each test and observed time series, p-values will 
be computed. In the framework of hypothesis test- 
ing under the null hypothesis, a p-value is the prob- 
ably of observing a value more extreme than the one 
actually observed. It allows to set a variability de- 
tection threshold independent of the number of mea- 
surements and of the noise level. 

The chosen tests are the test, tests on the asymmetry 
and kurtosis (these tests are related to three mo ments of 
the distribution). Abbe test ("von Newmannl ll941l) . a test 
to detect outliers, and a test on the presence of a slope. 

The stellar variability for certain variable types may re- 
main undetected in a wide band filter like the G-band be- 
cause of compensation effects in the spectral interval cov- 
ered by the G-band. The MBP and BBP bands will also 
be used. One method is to compute the covariant ma- 
trix between the different bands AIk,i — X]"=i(™fc,i ~ 
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™fc)(™/.i — rnl) / {<Jk.i<^i.i)7 where n is the number of 
measurements, the indices fc, I refer to the available pho- 
tometric bands and ml the mean Z-magnitude. The ma- 
trix diagonal terms are the values associated to the y^. If 
we diagonalise the matrix then we will get a signature of 
the physical nature of the variability. It seems an elegant 
way to combine the information of different filters and to 
enhance the variability detection capability. We will per- 
form some tests on the five Sloan filters of SDSS data and 
also on simulated data generated from synthetic spectra. 



6.2. Modelling the variability 

Once the variability is detected. Some further processing 
is require to characterise the variability. 

For strictly periodic behaviours, periodograms can be 
computed and the period of highest "power" can be used. 
There is no method which is superior to all the others for 
finding the period of a periodic signal, some are perform- 
ing better for certain types of behaviour, some others are 
more less prone to certain types of aliasing. The CPU 
consumption time may be also very different, though 
some optimisation are possible. We plan, as activities of 
the working group, to start benchmarking of these differ- 
ent methods. Simple method to recognise mono-periodic 
and multi-periodic phenomena should be also planed. 

Fourier series are adequate mathematical representation 
to model periodic variability (though often problematic 
for EA eclipsing binaries). However it appears diffi- 
cult to separate a set of well behaved light curves even 
after a Fourier deco mposition with analysis of residues 
(lEver & Blakel20Q4 . 

If the star is irregular, we have to rely on more general 
description: 

- The amplitude of variation: Several estimations are 
found in the literature. It is customary to give magnitudes 
at minimum and maximum light, however these quan- 
tities are not robust to outliers, and are meaningless if 
the amplitude of variation is small compared to the noise 
measurements. There are other estimations: signal vari- 
ance, inter-quartile (-decile) range, estimated amplitudes 
correct from noise measurements. 

- Time scales: The power spectrum can be used 
as a general description of the signal behaviour 
lEvans & Belokil^ (I200I. Calculation of variogram or 
auto-correlation algorithms maybe useful tools to provide 
time scales, when both the variable object and the time 
sampling are irregular 

These questions will be addressed. 



6.3. Classification 

We distinguish between two different strategies: a global 
approach to classification and an extraction of certain 
variable types. 



Mostly what has been done in the microlensing surveys 
is to extract some specific variability types. Criteria have 
been designed for know n variable typ e s, som e few ex- 
amples: RR Lyrae stars Alcock et al. (2000"), R Coro- 
nae Borealis Stars Alcock et al. (2001a), etc, from MA- 
CHO data, Cepheids l Udalski et al. 1999), eclipsing bi- 
naries (IWvrzvkowski et al..i.2004 ). etc, from OGLE data. 
On Hipparcos data, multivariate discriminant analysis 
has been applied to extract Slowly Pulsating B stars, 
and 7 Dor stars, cf. iWaelkens et aL C1998.) . .Aerts et all 

For the moment very few attempts of global clas- 
sification have been done, examples can be found 
m 'Poimanski & Macieiewski (2004) on ASAS data 
Ever & Blake (2004D on A SAS dat a, Belokurov et alJ 
(^03) on MACHO data, IWozniak et al.. (k2001l) . on 
ROTSE data, Marquette (private communication) on 
EROS data. 

We are therefore nowadays in a stage of method ex- 
ploration. Very promising methods like Self-Organising 
Maps (SOMs) are studied, see [Evans & Belokuroy 
( 2004) for a description of the method, cf. also Brett et al] 
(2004). Dr A. Naud (Torun Unversity) started a study on 
the global classification with SOMs of the Hipparcos pe- 
riodic variable star catalogue. This classification is using 
for each variable star its period, amplitude, colour V-I and 
the skewness of the photometric time series. The result 
is displayed in Figure |3] We observe that the classes of 
Mira, SRs, Cepheids and eclipsing binaries form isolated 
groups. The smaller amplitude variables are somewhat 
mixed. 

Several classification methods should be studied and 
compared: their efficiency, their capability to separate 
classes, their ability to detect new classes, their error lev- 
els, their robustness with respect to the sampling, their 
CPU consumption, and their level of required human in- 
teraction. 

Once some objects are recognised, further analysis can 
be done. Global treatment of eclipsing binaries hav e 
been done on OGLE data by Wvithe & Wilson (200 1^), 
IWvithe & Wilson (.2002) and are i nvestigated for the Gaia 
mission bv lPrsa & Zwitted(l2004 . 

Many tests can be run on existing surveys like ASAS, 
MACHO, EROS, OGLE. However we remind that the 
wealth of Gaia data is larger, the precision is higher and 
as for Hipparcos, the sampling is peculiar. Global classi- 
fication studies will have to be done with GDAAS. 

The global classification is not excluding procedures 
of specific tasks, selecting unusual and rare objects. 
As an example of an special selection procedure, 
we mention the detected HIPPARCOS planetary tran- 
sits of HP 2 09458, cf. Robichon & Arenou (2000) and 
ISod erhieTm ( 1999). It could have bee n singled o ut in Hip- 
parcos before the discovery o f.Charbonneau et al. (2000.). 
The HD 209458 HIPPARCOS photometric time series was 
not classified as variable, the eclipse depth is very shal- 
low with respect to the measurement noise, however with 
a well tailored statistical test ran on a subsample of Hip- 
parcos sources, HD 209458 is one of the best candidate 



6 



' »»••«•*•«••«••* 

rl^ ^.f^ « fli «>> # it ^1 «i « 



iiSiMte binaries 



*c* •••• ••• 

■•■«••» . « 
m ••»• ■ 
• 

• «*«• • •« 
•«■■•«■••■ 
•••>•••«•• 

i«« ' *»■•••«*»•■•• ' - ^ 
•«•*•'•«•»«'••• a*^ «•«•••«•■•••• >'•••«■«••'■••««•* 

• ••••• * *<•-•**••*••**••* «••••**»•*••••• 
••«*« #«**>•«■*«.•••••<••«•*>•«•»• 1 •••*•«•«»«•«•■• 

•«■■•• •«•••«•*•••••••••«•««• ••■•••'9a»*tt4*»H«»»*« 

••*•••**••*••• tia>*«***»*****«(*******«*«*»* 
••••«••«•••• ••••••9mK%4 •••••••••••a »< 

■» •••• •••••••••#l»«Mi»»MfM»»«f - ■■ ■ 

• ••#••**«••*«*•*«•«*••««••*•-•«*»• - '!' 

••*••••*•••#■•••*•«•••••••••••••••"> ' ' 

iSSaSSii*':- :• t' ••■••■»*««••■'':, 
••••*•**• ••*** •••***•*•*«••*- " 

*••+<•*•»• ***•«•*««•**«*««••••' 

■•«■*«■•••• «■••••••••«••«••••• 

••*••••*•••• •«•••••*■***«••**•* 

••• ••••» •«•••*•«•*«•«•••••• 

B« •*§» • ■»•■••••«••••••••• 
••••••••*•■«««««•• 
• ••• ••«• 

••«•* •*«• »•«•■•«••««•■•--; 



• ACYG 

• ACV 

BCEP 
BY 
CWA 
. CWB 
DCEP 
DCEPS 
DSCT 
EB 



••••••«•*•«• 

*•*•*•«•«■«• 
•••••••••••• 

1 *»•«••■«•*•••«••• 
>••••■••••«»•••• 



^ <••••••••••« 

. •>«•*•*•*••• 

• ••-••••«*••*•«•*»•«•••«•••«»•«•• 

• •§*••»♦»♦♦♦•••♦«•■»♦•»••♦««♦•♦♦»» 



A/B/Ci 



••»*« **•«««* 
•••*••«««••* 
. •■•*■••••••* 

«••*•••••••• 

«a»««*ft«tiii«« 

:ilRab' 



•••••• 

•«•• 



• «• ?*•«•••«••«•••«•••«••«••. ' «•••*••«« 

• •••* ".'#••«•••••»•*••«•••••••••*«* •»*•* 

a-«aa*aa*aaaa ^'•^taaaaaaaaaaataaaaaaaaaaaaav 

••••*••«•*••••«**•*•••••*••**••••«••*•*•*** 

•*aa'**a'*i^'^^9i# 

»»#aaa« ^ 



- -•9««aaaaf ' 
-.1 »«at - 
"-^a-aaiaaaaaaaaiiaaaaa* 



• EA 

• ELL 
« EVV 

E 

GCAS 

• M 

• RRC 

• RRAB 

• RS 

• SRB 

• SR 

• SARV 

• SRC 
SRA 

s SPB 

• SXARl 
S SXPHE 



*va'«**a«**«a*a^ 

:«••••**•*•*••• -!•«•*•«•«•««•• 
•■a*aa*aaaaaaa**#aaaaa 

' ••••••••••••• 

••••fIfnMug**? •»*••«••« 



----- 

li 

ii 



^»«»*aa*«a«»a*aaaa«a*l 

*•** « »«••«>*«•*•••••*••«••*•• 



SaMa 
•«aa 



••***M«**««M**M*M«**li*a*«*< 



Figure 3. Classification o/HIPPARCOS periodic variable 
star catalogue with Self-Organising Maps by Dr A. Naud. 



for having possible transits. 



7. CONCLUSION 



We have started to write algorithms for GDAAS to anal- 
yse the photometric time series. If the Gaia launch date 
of 20 11 seems far away, the tasks required to analyse the 
Gaia data for variability studies are enormous. If a gen- 
eral picture of the variability analysis is present, many 
problems are not solved yet. In order to gauge correctly 
the efficiency of proposed methods, of the analysis as a 
whole, there is the need to study thoroughly very techni- 
cal details. 

There are also broader questions: how and when to de- 
liver intermediate results on variable objects during the 
mission, how to plan Earth follow-up observations, how 
to combine Gaia data with coeval surveys, what tools to 
use for the visualisation of Gaia data and the quantities 
results of the variability analysis. 

If Gaia data will be an unprecedented data set, the vari- 
ability analysis also requires unprecedented ingenious 
methods, software developments and collaborative efforts 
of many researchers. 



ACKNOWLEDGMENTS 



I would like to express thanks to many members of the 
Variable Star Working Group, who are actively partici- 
pating to the tasks defined within the working group. Par- 
ticular thanks go to Dr. Dafydd Evans (co-taskleader) and 
Dr. Francesca Figueras, who are of continuous help and 
support. 

I would like also to thank for helpful comments and 
discussions Prof. M. Grenon, Prof. F. Mignard, Prof. 
D. Koch and Dr. F. Pont. 



REFERENCES 

Aerts, C, Eyer, L., & Kestens, E. 1998, A&A, 337, 790 

Alcock, C, AUsman, R., Alves, D. R., et al. 2000, ApJ, 
542, 257 

Alcock, C, Allsman, R. A., Alves, D. R., et al. 2001a, 
ApJ, 554, 298 

— . 2001b, ApJ, 562,337 

Belokurov, V., Evans, N. W., & Du, Y. L. 2003, MNRAS, 
341, 1373 

Belokurov, V. A. & Evans, N. W. 2003, MNRAS, 341, 
569 

Brett, D. R., West, R. G., & Wheatley, P J. 2004, MN- 
RAS, 353, 369 

Buzasi, D. 2002, in ASP Conf. Ser. 259: lAU Colloq. 
185: Radial and Nonradial Pulsationsn as Probes of 
Stellar Physics, 616-H- 

Charbonneau, D., Brown, T. M., Latham, D. W., & 
Mayor, M. 2000, ApJ, 529, L45 

ESA. 1997, ESA SP-1200 

Evans, N. W. & Belokurov, V. A. 2004, in this volume 

Eyer, L. & Blake, C. 2004, astro-ph/0406333 

Eyer, L. & Cuypers, J. 2000, in ASP Conf. Ser. 203: lAU 
Colloq. 176: The Impact of Large-Scale Surveys on 
Pulsating Star Research, 71-72 

Eyer, L. & Wozniak, P R. 2001, MNRAS, 327, 601 

Favata, F 1999, Baltic Astronomy, 8, 309 

Gott, J. R., Juric, M., Schlegel, D., et al. 2003, astro- 
ph/03 10571 

Ivezic, Z., Goldston, J., Finlator, K., et al. 2000, AJ, 120, 
963 

Ivezic, Z., Lupton, R. H., Anderson, S., et al. 2003, Mem- 
orie della Societa Astronomica Italiana, 74, 978 

Jordi, C. 2004, in this volume 

Pojmanski, G. 1997, Acta Astronomica, 47, 467 

Pojmanski, G. & Maciejewski, G. 2004, Acta Astronom- 
ica, 54, 153 

Pourbaix, D., Platais, I., Detournay, S., et al. 2003, A&A, 
399, 1167 

Prsa, A. & Zwitter, T. 2004, in this volume 
Robichon, N. & Arenou, F 2000, A&A, 355, 295 



7 



Schechter, P. 2005, ASP Conf. Ser. x: lAU Symp. 225 

Soderhjelm, S. 1999, Informational Bulletin on Variable 
Stars, 4816, 1 

Totani, T. & Panaitescu, A. 2002, ApJ, 576, 120 

Udalski, A. 2003, Acta Astronomica, 53, 291 

Udalski, A., Pietrzynski, G., Szymanski, M., et al. 2003, 
Acta Astronomica, 53, 133 

Udalski, A., Soszynski, I., Szymanski, M., et al. 1999, 
Acta Astronomica, 49, 437 

Udalski, A., Szymanski, M., Kubiak, M., et al. 2000, 
Acta Astronomica, 50, 307 

— . 2002, Acta Astronomica, 52, 217 

— . 1998, Acta Astronomica, 48, 147 

von Newmann, J. 1941, Annals of Math. Stat., 12, 367 

Waelkens, C, Aerts, C, Kestens, E., Grenon, M., & Eyer, 
L. 1998, A&A, 330, 215 

Walker, G., Matthews, J., Kuschnig, R., et al. 2003, 

PASP, 115, 1023 

Wozniak, P R., Akerlof, C, Amrose, S., et al. 2001, Bul- 
letin of the American Astronomical Society, 33, 1495 

Wozniak, P R., Udalski, A., Szymanski, M., et al. 2002, 
Acta Astronomica, 52, 129 

Wyithe, J. S. B. & Wilson, R. E. 2001, ApJ, 559, 260 

— . 2002, ApJ, 571, 293 

Wyrzykowski, L., Udalski, A., Kubiak, M., et al. 2004, 
Acta Astronomica, 54, 1 

Zebrun, K., Soszynski, I., Wozniak, P. R., et al. 2001, 
Acta Astronomica, 51, 317 

Zwitter, T. 2002, in ASP Conf. Ser. 279: Exotic Stars as 
Challenges to Evolution, 31-+ 

Zwitter, T. 2003, in ASP Conf. Ser. 298: GAIA Spec- 
troscopy: Science and Technology, 329-+ 



